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Based on standard perturbation theory, we present a full quantum derivation of the formula for 
the orbital magnetization in periodic systems. The derivation is generally valid for insulators with 
or without a Chern number, for metals at zero or finite temperatures, and at weak as well as strong 
magnetic fields. The formula is shown to be valid in the presence of electron-electron interaction, 
provided the one-electron energies and wave functions are calculated self-consistently within the 
framework of the exact current and spin density functional theory. 



Magnetism is one of the most important properties of 
materials. Both spin and orbital motion of electrons can 
contribute to the total magnetization. While the spin 
magnetization can already be calculated from first prin- 
ciples with high accuracy by state-of-art methods such as 
the spin density functional theory (SDFT), the study of 
orbital magnetization is still in a comparatively primitive 
stage. 

A first difficulty arises from the fact that there is 
still no theoretically well established formula for calcu- 
lating the orbital magnetization of a crystalline solid. 
The non-locality of the orbital magnetization operator 
M = —er x v is the major obstacle to obtaining a closed 
formula for an extended periodic system. Recently, Xiao 
et al. [1] and, independently, Thonhauser et al. [1, 0| 
obtained an orbital magnetization formula which avoids 
the non-locality problem and looks very promising for 
applications. However, up to date, there exists no gen- 
eral quantum mechanical derivation of this formula. The 
derivation presented in Ref. pj relies on the semi-classical 
wave-packet dynamics of Bloch electrons [H, 0, 0] , and its 
validity in the quantum context is not completely clear. 
On the other hand, the derivation presented in Ref. Q is 
quantum mechanical, but relies on the existence of local- 
ized Wannier functions, and cannot be easily generalized 
to metals or insulators with non-zero Chern number. In 
addition, both derivations are limited to non-interacting 
systems. The shortcomings of these approaches call for 
a full quantum mechanical and many-body theory of the 
orbital magnetization. 

A second difficulty is that a first principle calculation 
of the orbital magnetization (taking into account many- 
body effects) should be based on the spin current den- 
sity functional theory (SCDFT) 0, rather than the con- 
ventional SDFT. Unfortunately, SCDFT has been hin- 
dered so far by the lack of reliable expressions for the 
magnetization-dependent effective potentials. This may 
partly explain why the orbital moments of ferromagnetic 
transition metals such as Fe, Co, and Ni calculated in 
SCDFT were found to be significantly smaller than the 
experimentally determined values [ll[. How problematic 
these calculations are is well explained in the review ar- 
ticle by Richter The situation, however, has been 



rapidly changing in recent years. The advent of opti- 
mized effective potentials [la . [lij which treat exchange 
exactly and may systematically be improved for correla- 
tions opens new avenues to the study of magnetic mate- 
rials. Applications to atoms and molecules have already 
appeared in the literature, and applications to periodic 
systems are the obvious next step. 

Against this background, the present Letter serves a 
dual purpose. First, we present a general derivation 
of orbital magnetization in periodic systems based on 
the standard perturbation theory of quantum mechan- 
ics. The derivation clarifies the origin of the novel as- 
pects of the semi-classical derivation, such as the Berry 
phase correction to the density of states. It is gener- 
ally valid for metals and insulators with or without a 
Chern number, at zero or finite temperatures, in weak 
or strong magnetic field, and, of course, in the presence 
of spin-orbit coupling. Second, and most important, we 
combine this derivation with the exact current and spin 
density functional theory [f| 0], proving the validity of 
the magnetization formula for interacting systems. We 
believe that the magnetization formula, in combination 
with the recent advances the construction of optimized 
effective potential for SCDFT, will turn out to be a pow- 
erful practical tool for the study of systems that have 
long defied traditional ab-initio methods. 

We start from the standard thermodynamic definition 
of the orbital magnetization density: 
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where fl = E — TS — is the grand thermodynamic 
potential, V is the total volume of the system, and B 
is a magnetic field that only couples to the orbital mo- 
tion of electrons (but does not contribute to the Zeeman 
energy) [l^]. For convenience of derivation, we will first 
calculate the auxiliary quantity 
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where K = E - fiN. We have M = M + T (dS/dB) T ^ 
and, making use of the Maxwell relation (dS/dB) T = 
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(dM /dT) ^ B , we have the simple relation back to the 
orbital magnetization: 
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= M, 
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where /? = 1/kT. 

In principle, one can evaluate by employing the 
standard perturbation theory of quantum mechanics to 
calculate the energy correction due to a uniform mag- 
netic field. However, such an approach will again hit the 
difficulty of the non-locality of the orbital magnetization 
operator. To go around the difficulty, we apply an ex- 
ternal magnetic field that has an infinitely slow spatial 
variation 1 8, 91: 



B(r) = B cos(qy) z 



(4) 



The slow spatial variation of the field is controlled by the 
wave vector q, which will tend to zero at the end of the 
calculation. The correction to the energy density in this 
situation can be written, up to linear order in B, as 



5K(r) 



M ■ B(r) 



(5) 



We can then read out M, and take the limit of q — > 
for a uniform magnetic field. 



Non-interacting periodic systems — For clarity, we first 
carry out the perturbation calculation for non-interacting 
periodic systems. The single-particle Hamiltonian can 
be expanded as H = Hq + Vb, where Ho is the un- 
perturbed Hamiltonian, which yields the band disper- 
sion e„fc and the corresponding Bloch wave function 



ipnk(r) = exp (ik ■ r) u n k{r) , and Vb denotes the cou- 
pling to the external magnetic field: 



V B = - [t> • A(r) + A(r) ■ v] 



(6) 



where v is the velocity operator, and A(r) is the vector 
potential 



A(r) 



-B x, 



(7) 



which corresponds to the magnetic field discussed ear- 
lier. It is natural to define the grand-canonical ensemble 
energy density as K{r) = £ nfe f n kRe{^ k (r)Kip nk (r)}, 
where K = H — fiN, and f n k are the occupation number 
of the single-electron states of band index n and crystal- 
momentum k [19]. To first order in the perturbation, 
three kinds of terms arise from changes in the occupa- 
tion number, the operator K, and the wave function: 

5K{r) = Re^2 {sf nk ^ k K Q Tp nk + f n kipnk v B^nh 

+ fnk (lp* nk K 5lp nk + S^nkko^nkj } . (8) 

The orbital magnetization can be determined from the 
appropriate Fourier component of the energy density, i.e., 
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Ar5K(r) cos qy. 
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It is easy to verify that the first two terms of Eq. (JHJ 
does not contribute to M z , and only the contribution 
from the change of the wave functions remains. The first 
order perturbation to the wave function reads: 



e »(fe+9)-f \ Un , k+q ) (u n < k+q | v x (k) + v x (k + q) | u nk ) 



(10) 



where v(k) = dHo(k)/d(hk) is the velocity operator, and Ho(k) is defined from the unperturbed Hamiltonian by 
shifting the momentum operator with hk. The transformed Hamiltonian acts on the periodic functions, with u nk 
being its eigenfunctions and the band energy e nk its eigenvalues. Making use Eqs. ([8 UT0]) . we have: 
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Taking the long-wavelength limit q — ► 0, we obtain: 
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where f' nk = df(e nk )/de nk . The second term comes from the intra-band contribution with n = n'. Eq. lfT2]) can be 
further simplified with the help of the relations (u n > k | v x {k) \ u nk ) = (l/h)(e nk — e n > k ) (u n > k \ du nk /dk x ) for n ^ n', 
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and Im[. ..] = (1 / h)Im(du n k / dk y \e n (k) — Ho(k)\du n k/dk x ), where [. ..] denotes the expression inside the square 
bracket in Eq. lfT2]) , Combining these relations, and generalizing the result to the other components of M, we obtain 
finally: 
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The auxiliary and proper orbital magnetization be- 
come the same at zero temperature. In this case, the 
second term in the above expression Eq. I|13p vanishes 
because /' becomes a <5-function of (e n k — m)- The result 
is in perfect agreement with the semiclassical formula of 
zero temperature orbital magnetization of Xiao et al. [l[. 
For finite temperatures, we integrate Eq. J3]) and obtain: 
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(14) 

where m n (k) = (e/2ft)i(V fe u„ fc |[e n (fc) - H (k)} x 
VfcU„fc) is the orbital moment of state n, k and ft n (k) 
i (VfcU„fc | x | VfcU„fc) is the Berry curvature. The same 
expression was also obtained in Ref. [10]. 

Thus, all the previously known results are recovered 
by our fully quantum mechanical formalism. These 
results are valid not only for insulators with or without 
a Chern number, but also for metals at zero or finite 
temperatures. This implies that the semi-classical results 
are in general valid to linear order in the external fields. 
In hindsight, this should have been expected, because 
the semiclassical theory is designed to be exact in the 
limit of long length scales in the perturbation to the 
Hamiltonian. In our case, this length scale (through the 
vector potential) does diverge in the limit of vanishing 
magnetic field. 

Generalization to interacting systems — It is very de- 
sirable to generalize the above results to an interacting 
system. This can be done exactly within the framework 
of the current and spin density functional theory (CS- 
DFT) 0,0|. CSDFT is a generalization of the spin den- 
sity functional theory which includes the current density 
as an independent variable for the energy functional and 
thus provides direct access, via a variational principle of 
the Hohenberg-Kohn type, to the current density and the 
orbital magnetization of the thermodynamic equilibrium 
ensemble |Tl|. Following the formalism of Ref. 0], the 
many-body problem can be reduced to solving an effec- 
tive one-body Schrodinger equation: 



-L(-ihV + eA' a (r)f + V^r) 
2m 



with 



(15) 



K =V a + V H + V xca + — (Al- A'*) , (16) 
2m 

A' a =A a + A xc<7 . (17) 



Here V a and A a are the external scalar and vector po- 
tential, respectively, acting on the a component of the 
spin [a =| or j); Vr = e 2 J dr'n(r')/\r — r'\ is the 
Hartree potential, and V Kca and A Kca are the exchange- 
correlation (xc) scalar and vector potentials derived 
from the xc energy functional O xc [no-, j pa ] according to 
the formulas V xca - = 8Vl xc /8n a , eA XCIT = 5n xc /5j pa . 
The density, n,j(r), and the paramagnetic current den- 
sity, j P a( r ) are to be determined self-consistently from 
the solutions of the above equation according to the 
formulas n a (r) = £\ \^iv(r)\ 2 f{e ia ) and j pa {r) = 
(-a/MEitCMWfeM - V^ CT (r)^ CT (r)]/(e ia ). 
At finite temperature, the thermodynamic potential 
functional can be written as 0]: 

n = ~X> [l + e-^-/*)] 

ia 

}, (18) 

which is a functional of four fields: the densities n a and 
jpa , and the external potentials V a and A a - the last two 
enter the expression through the eigenvalues €i a . 

To calculate the orbital magnetization one needs to 
evaluate the variation SQ of the thermodynamic potential 
in response to a variation of the external magnetic field, 
which, in turn, is generated by a variation in the external 
vector potentials 5A a . In general, Sfl can be separated 
into two contributions: the primary one (Sfl\ n ,,j ptr ) arises 
directly from the variation of the vector potentials, keep- 
ing n a and j pa constant at their unperturbed values; the 
secondary one (<5S1| Aa ) might arise from the changes of 
n a and j pa at constant external potentials (these changes 
would affect ft via the modification of the effective po- 
tentials Vh, Vxco- and A xca ): 



sn = scii 
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It is easy to see that Sfl\ na j pcj contributes an orbital 
magnetization that is exactly given by Eq. lfl4]) , as if the 
system were a noninteracting system with eigenfunctions 
and eigenvalues determined by Eq. lfT5]) . This is because 
the variation of the external vector potential affects only 
the eigenvalues e ia in the one-body term of Eq. 1(18)1 . 
To evaluate 5Vt\y a ^A„ it is sufficient to observe that for 
given external potentials the thermodynamic potential Q, 
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is stationary against small changes of the density and 
the current about their equilibrium values: SVL/Sn a — 

sn/5j pa = o. 

Thus we have 

M\v.,A„ = °- (20) 

We then conclude that in the context of the CSDFT, 
we can treat the system as an effective one-body system, 
and use Eq. lfl4|) to calculate the orbital magnetization, 
albeit using the dispersion and wave-functions derived 
from Eq. lfl5|l . 

We stress that this conclusion would not hold true 
if the one-electron orbitals and their energies were 
calculated within the framework of the ordinary density 
functional theory [l5| (as opposed to CSDFT). In such 
a formulation the xc energy functional would be a 
functional of density and magnetic field: Q xc [n a , B]. 
Then the formula for the orbital magnetization would 
include a functional derivative of f2 xc with respect to B, 
at variance to the one-electron formula of Eq. lfi"4"ll . 

Finite fields — Our formula can be applied rigorously 
to finite magnetic fields if these fields are rational in the 
sense that fluxes through the faces of a unit cell are frac- 
tional multiples of the flux quantum h/e. In this case, one 
can define Bloch like eigen-states with respect to mag- 
netic translations, which are ordinary translations on the 
crystal combined with gauge transformations. Orbital 
magnetization at such a field can then be calculated per- 
turbatively by adding a small change SB to this field. 
Both the semiclassical theory and quantum perturbation 
with respect to SB give the same expression for the or- 
bital magnetization, provided we use the magnetic Bloch 
wavefunctions as a basis. Like wise, the CSDFT can also 
be formulated for systems with periodic boundary condi- 



tions with respect to magnetic translations, so the justifi- 
cation of our results for interacting systems is straightfor- 
ward Q- Indeed, current density functional theory has 
been used to study the formation of an electron crystal 
(Wigner crystal) at very high magnetic field |lGl |. 

The situation for irrational fields is a bit tricky if 
one insists on rigorous results (l7l |. One may consider 
an irrational field as a limit of sequence of rational 
fields. This is possible if the magnetization depends 
on the field continuously, which is expected to be the 
case when the temperature is finite. Indeed, when the 
fast de Haas-van Alphen oscillations are smeared out, 
the average magnetization changes continuously with 
the Fermi energy [131 - For a fixed Fermi energy, the 
average magnetization also changes continuously with 
the magnetic field. Therefore, we expect that the average 
magnetization is a continuous function of magnetic field 
at a fixed density of electrons. 

In summary, we have presented a full quantum deriva- 
tion of the orbital magnetization formula. The derivation 
is generally valid for insulators with or without a Chern 
number, for metals at zero or finite temperatures, and at 
weak as well as strong magnetic fields. We also find that 
the resulting formula is directly applicable to interact- 
ing systems provided one uses one-electron energies and 
wave functions obtained from the self-consistent solution 
of the Kohn-Sham equation of current and spin-density 
functional theory. 
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